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We study the glueball properties at finite temperature below T c using SU(3) anisotropic quenched lattice QCD 
with j3 = 6.25, the renormalized anisotropy 7 = a s /a t = 4 and 20 3 x N t (N t = 35,36,37,38,40,43,45,50,72). 
From the temporal correlation analysis with the smearing method, we observe the mass reduction of about 20 % 
for the lowest ++ glueball as m G (T) = 1.25 ± 0.1 GeV for 0.8T C < T < T c in comparison with m G ^ 1.5 ~ 1.7 
GeV at T ~ 0. 



1. INTRODUCTION 

Finite temperature QCD, including the quark 
gluon plasma(QGP) physics, is one of the most 
interesting subjects. Even in the hadronic phase, 
a lot of effective models suggest the changes in 
the hadronic properties. However, lattice QCD 
studies of thermal hadron properties have been 
inadequate until quite recent. This is due to the 
difficulty in measuring hadronic two-point corre- 
lators at finite temperature. At high tempera- 
ture, due to the shrink in the temporal exten- 
tion and the consequent decrease in the number 
of data, the mass measurement from the tempo- 
ral correlations is difficult. Hence, finite temper- 
ature mass shifts have been studied through the 
spatial correlations, which, however, is afflicted 
with the mixture of large Matsubara frequencies 
jlj. Recently, the use of anisotropic lattice is es- 
tablished |Q , which has the finer temporal lattice 
spacing at than the spatial one a s . In this way, 
available data increases, and the accurate mass 
measurement directly from the temporal correla- 
tions becomes possible |3|-|| ■ Here, we report the 
thermal properties of the ++ glueball based on 
the quenched SU(3) anisotropic lattice QCD. 



"The lattice calculations have been performed on NEC- 
SX5 at Osaka University. 



2. SMEARING METHOD 

We consider the glueball correlator as G(t) = 
(O(f)O(O)), O(f) = ^(0(£,f)-(6)) where 

X 

6(x,t) = ReTr(P 12 (x,t)+P 23 (x,t)+P 31 (x,t)) de- 
notes an interpolating field for the ++ glueball. 
Pij(x,t) € SU(3) denotes the plaquette operator 
in the i-j-plain. We apply the spectral represen- 
tation to O(f) to have G(V)/G(0) = ^C^e"^*, 
where C n = \ (n\O\0) | 2 / £ | (fc|O|0) | 2 with E n 
being the energy of the n-th excited state \n). 
Here, |0) denotes the vacuum and |1) denotes 
the ground state glueball. C n is non-negative 
with ^C„ = 1. It is known that G{t) re- 
ceives large contributions from the excited states 
and the ground state contribution is small. As 
a consequence, the extracted mass always be- 
haves much heaver than the ground-state mass. 
The problem originates from the fact that the 
"size" of the ordinary plaquette operator is 0(a s ), 
which is much smaller than the physical size R 
of the glueball. Hence, the problem can be re- 
solved by generating physically extended glue- 
ball operators, which is achieved by the smear- 
ing method Q. This extended operator is ref- 
crcd to as the smeared operator, which is ob- 
tained by replacing the ordinary link variables 
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Ui(s) in the plaquette operator Pij(s) by the 
fat link variables J7,(s) S SU(3). L/j(s) is de- 
fined so as to maximize RcTr (u] (s)Vi (s)^ , where 

Vi(a) = aUi(s) + U ±j (a)U i (a±j)ui j {a + i) 

with U-^(s) = 11^(3 — p,). a is a real parameter. 
The summation involves only the spatial direc- 
tions to avoid the temporal nonlocality. Note that 
Ui(s) holds the same gauge transformation prop- 
erties as Ui(s). We refer to the fat link defined in 
this way as the first fat link U^(a). The n-th fat 
link U- n \s) is defined iteratively by U- n \s) = 
I7.-™ _1) (s) starting from C/f (s) = U l (s). The 
plaquette operator constructed with (s) is 
refered to as the n-th smeared operator. We next 
consider the size of the n-th smeared operator. 
By using the linearization and the continuum ap- 
proximation, we obtain the diffution equation as 

^-<f> i (n;Z)=DA4 >i (n;x), D=-^-. (1) 
on a + 4 

4>i(n;x) describes the distribution of the gluon 
field Ai(x,t) in the n-th smeared plaquette. In 
the n-th smeared plaquette located at the origin 
x ~ Q, the gauge field Ai(x, t) is distributed in the 
Gaussian form as e~ s2 /{iDn) / (A-KDnf 2 . Hence 
the size of the operator is estimated as 



R=^) = V6D^ = aJ^j. (2) 

V a + 4 

The original aim of the smearing method is the 
accurate mass measurement by maximizing the 
ground-state overlap. However, it can be also 
used to give a rough estimate of the physical glue- 
ball size. In fact, once the maximum overlap is 
achieved with some n and a, the glueball size can 
be estimated with Eq.(||). 

3. SU(3) LATTICE QCD RESULT 

We use the SU(3) anisotropic lattice plaquette 
action as S G = -^-i RcTr(l - Pij(a)) + 

C s,i<j<3 

j^Co ReTr(l — _P;4(s)) with the plaquette op- 

c s,i<3 

erator P^ v (s) in the /i-^-plane. The lattice pa- 
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Figure 1 . The 0++ glueball correlator G(t) /G(0) 
for N BmeaI = 40 at T = 245 MeV {N t = 37). The 
statistical errors are hidden within the symbols. 
The solid line denotes the best single hyperbolic 
cosine fit performed in the interval indicated by 
vertical dotted lines. 



rameter and the bare anisotropy parameter are 
fixed as f3 — 6.25, £0 = 3.2552, respectively 
These parameters reproduce the spatial lattice 
spacing as a" 1 = 2.272(16) GeV, and the tem- 
poral one as a^ 1 = 9.088(64) GeV. The rcnor- 
malized anisotropy is 7 = a s /a t = 4 Here, 
the scale unit is determined by reproducing the 
string tension as \fa = 427 MeV from the on- 
axis data of the static inter-quark potential. The 
pseudo-heat-bath algorithm is used to update the 
gauge field configurations on the lattice of the 
sizes 20 3 x N t with various N t listed in Tabic |l|. 
For each temperature, we pick up gauge field con- 
figurations every 100 sweeps for measurements, 
after skipping more than 20,000 sweeps of the 
thcrmalization. The numbers of gauge configu- 
rations used in our calculations are summarized 
in Table |l|. From the analysis of the Polyakov 
loop, we find T c ~ 260 MeV. To enhance the 
ground-state contribution, we adopt the smear- 
ing method with a — 2.1. 

In Fig.|, the 0++ glueball correlator G(t)/G(0) 
for iVsmear = 40 at T = 245 MeV is shown. 
The statistical errors are estimated with the jack- 
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Table 1 

The lattice QCD result for the lowest ++ glueball mass at finite temperature. The temporal lattice 
points Nt, the corresponding temperature T, the lowest ++ glueball mass rriQ(T), the maximal value 
of the ground-state overlap C max , uncorrelated x 2 /Nbf, the best smearing number iV^f ar , the smearing 
window N^t? w > the number of gauge configurations V con f and the glueball size R ~ (V^) are listed. 



N t 


T [MeV] 


m G [GeV] 


^max 


x'/Nuf 


best 
smear 


Afwindow 
smear 


^Vconf 


R [fin] 


72 


125 


1.54(1) 


0.981(6) 


0.23 


41 


32 - 


51 


1284 


0.49 ~ 0.61 


50 


182 


1.43(2) 


0.967(8) 


0.02 


45 


36 - 


55 


1051 


0.52 - 0.64 


45 


202 


1.48(2) 


0.971(8) 


0.08 


37 


30- 


46 


1382 


0.47-0.58 


43 


210 


1.30(2) 


0.945(8) 


0.35 


40 


32 - 


52 


1197 


0.49 - 0.62 


40 


227 


1.29(2) 


0.935(6) 


1.22 


40 


32 - 


51 


2021 


0.49 - 0.61 


38 


239 


1.36(2) 


0.935(7) 


1.00 


38 


30- 


50 


2030 


0.47-0.61 


37 


245 


1.18(3) 


0.889(12) 


0.91 


43 


35 - 


53 


2150 


0.51 - 0.63 


3G 


252 


1.38(2) 


0.948(8) 


0.09 


36 


29 - 


46 


1744 


0.46 - 0.58 


35 


259 


1.26(2) 


0.931(7) 


1.15 


40 


31 - 


51 


2244 


0.48 - 0.61 



knife analysis. The solid line denotes the best 
single hyperbolic cosine fitting as G(t)/G(Q) = 
C ( e -m G ta t + e -m G (N t -t)a t ^ per f orme d i n the in- 
terval indicated by vertical dotted lines. In 
the most suitable smearing number A/^|^ r , the 
ground-state overlap C is maximized and the 
mass tog is minimized, which indicates the 
achievement of the ground-state enhancement. In 
practical calculations, the maximum overlap and 
the mass minimization are achieved at almost 
the same iV smear , and both of these two condi- 
tions would work as an indication of the max- 
imal ground-state enhancement. Here, we take 
the maximum ground-state overlap condition. 

In Table g, we summarize the SU(3) lattice 
QCD result of the lowest ++ glueball mass 
tog(T) at various temperatures. We observe that 
about 20 % mass reduction of the lowest ++ 
glueball near T c as to g (T) = 1.25 ± 0.1 GeV for 
0.8T C < T < T c in comparison with m,Q(T — 0) — 
1.5-1.7 GeV |0,§. To estimate the glueball size, 
we search iViiS!,* which achieves the maximum 

smear 

overlap as C = C max . Considering about 0.5 
% errors of C, we define the "smearing window" 
N S m"ar OW where C > 0.995C max is satisfied. From 
Eq.(g) with N™f° w , we estimate the glueball size 
as R — 0.5 — 0.6 fm both at low and high temper- 
atures. In Table |l|, we summarize the lowest ++ 
glueball mass toq(T), the ground-state overlap 
C max , uncorrelated xV^df, N^, A&° w , 
the number of gauge configurations N con { and the 



glueball size R. 

4. SUMMARY 

We have studied the ++ glueball properties 
at finite temperature using SU(3) anisotropic 
quenched lattice QCD with more than 1,000 
gauge configurations at each temperature. From 
the temporal correlation analysis with the smear- 
ing method, we have observed about 20 % mass 
reduction of the lowest ++ glueball as ma(T) = 
1.25 ± 0.1 GeV for 0.8T C < T < T c , while no sig- 
nificant change has been observed for the glueball 
size as R — 0.5 — 0.6 fm. 
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